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Abstract 

To predict the heat diffusion in a given region over time, it is often 
necessary to find the numerical solution for heat equation. With the 
techniques of discrete differential calculus, we propose two uncondi- 
tional stable numerical schemes for simulation heat equation on space 
manifold and time. The analysis of their stability and error is accom- 
plished by the use of maximum principle. 
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1 Introduction 

To investigate the predictions of heat equation it is often necessary to ap- 
proximate its solution numerically [HE]- Discrete exterior calculus (DEC) 
constitutes a discrete realization of the exterior differential forms structure, 
and therefore, the right framework in which to develop a discretization of 
exterior differential system |3l-|10j. The operators in this system such as 
Hodge star, exterior derivative, and Laplace operator can also be naturally 
discretized by DEC. The techniques of DEC can be used to construct a con- 
ditional stable scheme for heat equation [11] . The stable condition is a very 
severe restriction, and implies that very many time steps will be necessary 
to follow the solution over a reasonably large time interval. 
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In this paper we shall propose two unconditional stable schemes for the 
numerical solution of heat equation in space manifold and the time, namely 
implicit and semi-implicit DEC schemes. The analysis of their error is ac- 
complished by the use of maximum principle. The methods proposed here 
can be extended to problems with general boundary condition as explicit 
scheme, then to general linear and nonlinear parabolic equations. 

2 Implicit DEC scheme for heat equation 

For a function ip{x,y,z,t) of three spatial variables {x,y,z) and the time 
variable t, the heat equation is 

with boundary conditions where ip is the temperature, p is the material 
density, c is the material specific heat, k is the thermal conductivity, Q is 
the internal heat source density. 

The DEC method can approximate Laplace operator as 

A *'^(f * +<f * d, 

where d is the matrix of discrete exterior derivative, * is the matrix of 
discrete Hodge star operator. The explicit DEC scheme for heat equation 
uses the forward time difference for the temporal derivative of Eq.(l), which 
is conditional stable [llj . If we need to refine the space mesh to improve the 
accuracy of the solution the amount of work involved increases very rapidly, 
since we shall also reduce the length of time step. Now, we show how the 
use of a backward time difference gives a difference scheme which avoids this 
restriction. 

Given 2D or 3D space manifold with a smooth boundary, a simplicial 
mesh is a tessellation of it by triangular or tetrahedra, satisfies the condition 
that any two of them may intersect along a common face, edge or vertice. 
For some situations, a source having azimuthal symmetry about its axis is 
considered. In this case, we only need to consider 2D triangular discrete 
manifold as the space. The 3D case can also be done in the same way. 

At first, we define some values on mesh. Take Fig.l as an example for 
2D mesh, in which 0, A,..., F are vertices, 1, 2,. ..,6 are the circumcenters of 
triangles, a, b,...,f are the circumcenters of edges. Denote kj as the length 
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of line segment (i, j) and Aijk as the area of triangle (i, j, k). 



Define 
and 



Figure. 1 

^12 := hf + hf, '■= ha + ha, •■•) '■= he + he. 



-Pl23456 ^01/ + ^02/ + • • • + ^Ole- 

The Laplace operator acting on tp at vertice approximated by DEC is 

A^/'o ~ -s- — ( r^lV'vi - V'o) + 7^(V'B - ^o) + 7^(^c - V'o) 

A23456 V'AO 'BO 'CO 

+7^(V'D - V'o) + -T-i.'^E - i'o) + 7^(V'F - V'o) I • 

I'DO t_EO tFO / 

If we use the backward time difference, we obtain the implicit DEC scheme 
for Eq.(l) as follows: 

C/J-ri23456 'BO 
'CO 'DO 'BO 



(2) 

Now, we use maximum principle to prove the unconditional stability for 
scheme (2). Given a perturbation e" on each the relation between e" 
and e"^^ can be induced from scheme (2) as follows: 

n+l _ At / /23 ^34 . ^45 , ^56 , ^16 . ^12 \ ^„+l 



° ° -F'l23456 V^O 'bo 'go 'bO 'bO 'bO / ° 

At / ^23 n+l . hA «+l . ^45 «+l . ^56 «+l 
'T— 1: A -f- En -f- Er<. ~T 1 £; 



A23456 V^Ao ^ Ibq ^ Icq 'bo ^ 



h& n+l I '12 n+l 

+7 ^B 
tBO 'BO 
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If we introduce the maximum perturbation at n time step by writing 

|e"| = Maxiey|e"|, 



we can obtain 



I n+ii ^ A23456 V^Ao Ibo Icq I do Ieo Ifo 

\£q I S 



At ^ l23_ _|_ h4_ _|_ ^ h6_ _|_ _|_ l^n+1 1 



-P123456 V^AO ^BO ^CO ^DO ^BO ^FO 



(3) 



Inequahty (3) holds for all vertices in grid, therefore 



That is to say scheme (2) is unconditional stability. Longer time step can 
be used than in explicit DEC scheme. 

By the definition of truncation error, the solution ijj of Eq.(l) satisfies 
the same relation as scheme (2) except for an additional term 0((At)^ + 1^) 
on the right hand side. Thus the error X" = ipf — ipf is determined from 
the relation 

-^n+l ^ j^n _ (^i^ + i3£_^j45__^j56__^_fl6__^ _fl2_\ 

° -P123456 V'ao Ibo ^co ^do Ieo Ifo J ° 



A23456 V^Ao Ibo ^ Ico ^ Ido ^ 



Define 

It follows that 
and hence that 
so that, since X° = 0, 



< |X"| + 0((At)2 + Z2)^ 

|j^n+l| < |X"|+0((At)2 + /2)^ 

< 0{{Atf + f), 



and this tends to along the refinement path under the assumed hypotheses. 
So far we have assumed that numerical errors arise from the truncation 
errors of the finite difference approximations, but that the boundary values 
are used exactly. 
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3 Semi-implicit DEC scheme for heat equation 



Implicit DEC scheme is not so easy to use as the exphcit DEC scheme. We 
must solve a system of equations to give the values simultaneously, which 
will cost a lot of time. Now we propose a semi-implicit scheme, which 
combines the virtues of explicit scheme (direct calculation) and implicit 
scheme (unconditional stability) . This scheme is that the calculated element 
using next point of time whereas the other elements using current point of 
time. Hence the semi-implicit scheme of Eq.(l) is 

tCO , 'do 'BO 



It can be rewritten as 



^"123456 V 'AO 'bo 'CO 



(4) 

Now, we consider the stability and convergence for scheme (4). The 
relation between perturbation and e"^^ can be induced from scheme (4) 
as follows: 

+ ( b^^n i34 ^„ ^ h5_^n i56_^„ _^ _fl6 ^„ ^ hi 

_ ° CpPi23456 \Iao ^ IbO ^ Ico ^ IpO ^ IeO ^ IpO ^ 

° 1-1- f bl. -J- _|_ j45_ he_ h6_ ^ h2 

cpPi23456 V^Ao Ibo Icq Ido Ieo Ifo, 

(5) 

Since the coefficients are positive, we can omit the modulus signs in the 
Eq.(5) to give inequality 

, kAt ( hz 

|g.n+l| ^ CpPi23456 \Iau 

^ \ -\- fcAt / ^23 _|_ ^34 _|_ Us _|_ ha _|_ Wa _|_ li2 ^ ' ' (6) 

C/9P123456 \Iao 

< \e'^\. 

Inequality (6) holds for all vertices in mesh, therefore 
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That is to say scheme(5) is unconditional stability. 

The relation of error can derived form scheme (4) as follows: 

]^ I fcAt ( '23 I '34 I '45 I 'sB I I '12 A 

cp^'123456 V 'AO 'sO 'CO 'dO 'eO 'fO / 
+0((At)2 + Z2). 

Because of the nonnegative coefficients, it follows that 

< |X"| + 0((At)2 + i2)^ 

and hence that 

|^n+l| < |X"|+0((At)2 + i2)^ 

SO that, since X° = 0, 

<0((At)2+z2)^ 

and this tends to along the refinement path under the assumed hypotheses. 
That is to say scheme (5) is convergence. 

In scheme (2) and (4), the derivative is approximated by first order dif- 
ference. Equivalently, is approximated by linear interpolation functions. 
Consulting the definition about accuracy of finite volume method, we can 
also say that scheme (2) and (4) have first order temporal and spacial accu- 
racy. 

4 Implementation of semi-implicit DEC scheme 

The implementation of semi-implicit scheme (4) for heat equation is the 
same as the explicit scheme consisting of the following steps: 

1. Set the simulation parameters. These are the dimensions of the com- 
putational mesh and the size of the time step, etc.; 

2. Initialize the mesh indexes. 

3. Assign source. 

4. Assign boundary conditions. 

5. Compute the value of all spatial nodes and temporarily store the result 
in the circular buffer for further computation. 

6. Visualize the currently computed grid of spatial nodes. 
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7. Repeat the process from the step 3, until reach the desired total num- 
ber of iterations. 



In the following examples, the parameters [sp = k = c= l and the heat 
source is y^t/SOO. The examples in Fig. 2 and Fig. 3 show that the diffusion 
of heat on surfaces of dragon and torus. The temperature -0 > 1 in red 
domain, < ^|J < 1 in yellow domain, and ip = in blue domain. 



Figure. 2 




Figure. 3 
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